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I. INTRODUCTION 



The detection of Higgs particles is one of the most important objectives of the Large 
Hadron Colhder (LHC). Charged Higgs bosons are predicted in extended versions of the 
Standard model (SM), like two-Higgs-doublet models (2HDM) and the Minimal Supersym- 
metric Standard Model (MSSM). A discovery of such an additional Higgs boson will imme- 
diately indicate physics beyond the SM; there is, hence, strong theoretical and experimental 
activity to provide the basis for its accurate exploration. 

At hadron colliders, the charged Higgs boson could appear as the decay product of 
primarily produced top quarks if the mass of iJ^ is smaller than irtt — nib. For heavier H^, 
other mechanism for H"^ production have been investigated: single Higgs-boson production 
associated with heavy quarks, like gh H't and qb —>■ q'hH~ |Q , and pair production of 
if^ through tree-level qq annihilation and via the loop- induced gluon-fusion mechanism . 
Moreover, single //^production in association with a W boson, via tree-level hh annihilation 
and one- loop gg fusion has been proposed and analyzed 0. Detailed studies show that 
these production mechanisms at the LHC can help to explore the parameter space, even 
beyond ~ ITeV and down to at least tan (3 ~ 3. 

Since U) annihilation is the main source of the hadronic W^H^ production cross section, 
it is necessary to calculate and implement also the loop contributions to hb W^H^ 
for more accurate theoretical predictions. Recently, the calculation of the 0{aew'fn't{h)l^w) 
and 0{aewmt{h)l^w) supersymmetric electroweak(EW) one-loop corrections were presented 
in in the frame of the MSSM, which can give rise to a 10-15% reduction of the lowest-order 
result. In this paper, we deal with the one-loop QCD corrections to hh W^H^ . 

The arrangement of this paper is as follows. Section II contains the analytic results, and 
in Section HI we present numerical examples and discuss the implications of our results The 
lengthy expressions of the form factors are collected in the Appendix. 

II. ANALYTIC EXPRESSIONS 
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A. Virtual corrections 



The FejTiman diagrams for charged Higgs-boson production via the parton process 
b{pi)b{p2) H'{ki)W^{k2), including the QCD corrections, are shown in Fig. 1. The 
diagrams are created by use of FeynArts [0] and are handled with the help of FeynCalc P] . 
We keep the finite 6-quark mass throughout the calculation, in order to control the coUinear 
divergences. 

As usual, we define the Mandelstam variables as 

s = {pi+P2f = {ki + k2)\ 
i = {pi - k^f = {p2 - /ca)', 

u={pi- k2f = {P2 - k,f. (1) 

Taking into account the 0{as) corrections, the renormalized amplitude for hb — > W^H~ 
can be written in the following way, 

M^'^ and M^^ are the s- and t-channel tree- level diagrams corresponding to Fig. 1(a) and 
Fig. 1(b), which are given by 



cosa sin(Q; — /?) sin a cos(q; — /?) 



s — my s — mj^ 



(Ms + Mg + Mq + M: 



10 J 



= ^ l^b sin /3(2M9 + m^Mi + M3) - htirit cos PM2] , (3) 

V2(t - mt) 



where 



, grrib gmt 

y/2mw cos p Vimwsmp 

denote the Yukawa couplings of the bottom and top quarks. Mj are reduced standard matrix 
elements, which are defined by 
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M2=V{P2) Mk2)PLuipi), 
M3=V{P2) M Mk2)PRu{pi), 

M^ = v{p2) M Mk2)PMpi), 
M5 = v{p2)Pru{pi)pi ■ e{k2), 
Mq = v{p2)Plu{pi)pi ■ e{k2), 
M7 = v{p2) hPRu{pi)Pi ■ eih), 
Ms = v{p2) hPLuipijPi ■ e{k2), 
Mg = v{p2)Pru{pi)p2 ■ e{k2), 
Mio = v{p2)Plu{pi)p2 ■ £{k2), 
Mu = v{p2) hPRu{pi)P2 ■ e{k2), 

M12 = v{p2) MPlu{pi)p2 ■ £{k2), (5) 

with the projectors Pl,r = (1 T 75)/2. 

The terms 6M in describe the virtual contributions from the 2- and 3-point functions 
together with their counterterms, and 6M^°^ denotes the contribution from the irreducible 
4-point function. The vertex and self-energy corrections to the tree-level process are included 
in 6M^'^ , which are given by 

c- ,",1/, C^W :r:r ,X 9^fer COS Q! SlU ( O " /5) SlU O COS ( O! — /3) , , , ^ ,^ ,^ ,^ , 

6M'''^'\H, h) = ^ . ^ 2 + Me + Mg + Mio 

V2 s — mjj s — mf^ 

,6mh 1 1 



mi, z, z, 



V2(s — mj^) nib Z z 



SM^iit) = ^9 ^j^^ sm(3{2Mg + mfeMi + M3) - htmt cos/JMa] 



V2{i 



m 



5M^W = J [hsmP{2Mg + m^M, + M,){2mt6mt + {m\ - t)Z\) 
y z[t — rrif ) 

-ht cos l3M2{5mt{ml + i) + rritiml - i)Zi)] + 5M^^*\ 
5M^^W = A h sin/?(2M9 + m^M, + M,){^ + \zi + \zl) 
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-htmt cos/3M2(^ + \zl + ^Z^)] + 5M^^W. (6) 

Therein, 5M^i(*)(if, /i, A), 5M^i(*), (5M"^W, (5M^2(i) ^present the unrenormahzed one-loop 
corrections arising, respectively, from the bbH{h,A) vertex diagrams in Fig. 1(c), the tbH~ 
vertex diagram Fig. 1(d), the btW^ vertex diagram Fig. 1(e), and the top-quark self-energy 
diagram Fig. 1(f); 6M^°^ corresponds to the box diagram Fig. 1(g). All the 6M^'^ and 
SM^°^ can be written in the form 

5M = CfY.UM,, Cf = \, (7) 

with the reduced matrix elements Mj in (^ and form factors /«, which are given explicitly 
in the Appendix. 



Here we present the analytical results in the MS scheme; accordingly, the quark masses 
in (^) are running masses rrib^tifJ')- The results in the OS scheme can be easily obtained 
by using the pole masses and replacing the corresponding field- and mass-renormalization 



constants. Actually in the MS scheme, we must add the finite part of the wave-function 
renormalization constants for the b quarks, according to the LSZ prescription. 

From the self-energy diagram in Fig. 1(h) for the b, t quarks we get the explicit expressions 
of the renormalization constants, valid for both the t and b quark, 

^ = -^3C,A, 
m Att 

cx 

Zl = Zr = —-r-CpA, with 
47r 

A = ^ -7i5 + log(47r), e = 4-D, (8) 

in dimensional regularization. 

Squaring the amplitude and performing the spin summations, 

E l^-r = E 1^^^ + ^o*^r + 2ReE'5M(Mi^) + M^^V , (9) 
yield the partonic cross section for the process bb W^H~ with virtual corrections. 



with 
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i± = + ' ± - (^^ + mH-y){s - {mw - mn-Y) ■ (11) 

The partonic cross section d^l^^y^_^+ioT the process hb W~H^ is identical to that of 
bb W+H-. 



B. Real-gluon contributions 

The virtual corrections involve an infrared singularity from the massless gluon. For 
our purpose, we can use a small gluon mass A for regularization, which allows to identify 
the infrared divergence as a log(A) term. This infrared divergence is cancelled by adding 
the corresponding real-gluon-radiation corrections, displayed in Fig. l(i-m). For technical 
reasons, it is convenient to perform the phase-space integration over a soft- and hard-gluon 
part separately. The soft-gluon contribution to the cross section at the partonic level is 
proportional to the tree-level cross section a^-^^^j^^^-^ for hb H^W^ , 

'soft _ - o_ r 

Ss = -C^-^{log — ^ + log + log log — ^ + ^ + o log (-:^)}' (12) 

ZTT a S S a O z s 



where AE is the energy cutoff for soft gluons. After adding the cross section (|ToD with the 
virtual corrections, the sum is independent on A, and the log^ term cancels. 

For the hard-gluon part, we use the Monte Carlo packages BASES to perform the 
phase space integration, with the cutoff that the energy of the gluon is greater than AE. We 
do not give detailed expressions here. Numerically it was checked that the sum 0"*^°^* -|- a^'^'^'^ 
is independent of the cutoff AE. 

At 0{as), there are also initial-gluon contributions from gb bW^H^, the Feynman 
diagrams can be obtained from Fig. l(i-m) by treating gluon and b quark as incoming 
partons. For the processes, we should caution how to subtract the on-shell top quark and/or 
charged Higgs boson contributions. Following the methods of Ref. |T^, the partonic cross 



section after subtracting the on-shell contributions of subprocess gb bW^H^ can be 
written as: 

^lbZlw±H^ = ^gb^bw^H^ - ^gb->tH-Br{t bW+) 



—a 





gb^tW 



^Br(t bH+). (13) 



The cross section of gb bW^H^ can be obtained similarly. 

C. Hadronic cross section and subtracting double counting 

In the approach described above we have considered the real-gluon corrections, which 
give rise to a term involving log(s/m^) from the region where the b quark splits into a 
collinear 6-quark-gluon pair and the gluon into a collinear b pair. This logarithm is already 
contained in the heavy-quark distribution function, hence it has to be removed. This can be 
done, following [0, by subtracting the tree-level process bb H^W^ convoluted with one 
heavy-quark distribution function given by the perturbative solution to DGLAP equation, 

zn rrijj Jx y y 

= ^T— log(-i) / —Pqg{-)9{y,f^f), (14) 
/TT Jx y y 

where 

P,,{z) = \[z' + {l-zf]. (15) 
In this way, the total cross section at 0{as) can be expressed as 

= b®a^l®b- b(^bb) ® &u^h±w^ ® 6 - 6 ® '^56_.i^±w'T ® kbb) 
+b d-'^bgZlH±w^ 0g-b0 d-I^^H±w^ ® kbg) 

+b a^^^j^^^^^ ®g-b® a^Vz/iiyT ® kbg) , (16) 

where a^i is the infrared-finite parton cross section which is given by the sum of a™^ + d^'^^^ + 
^hard_ (p^G]), A® (J ® B represents the cross section of the subprocess a is convoluted 

with the parton distribution functions (PDF) A and B, 

A®a®B= / dz — a{z^s) , zq = , (17) 

J dz Y ^ 

where ^/s is the overall CM energy of the pp system, and the parton luminosity dL/dz is 
defined as 

— = 2zJ^^-A{x,^f)B{-,^f). (18) 
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III. NUMERICAL RESULTS 



We now present a numerical discussion of the QCD corrections to W^H"^ associated 
production at the LHC The SM input parameters in our calculations were taken to be 
aew{mz) = 1/128.8, as{mz) = 0.118, mw = 80.41GeV and mz = 91.1867GeV, and the 
pole masses of top and bottom quarks are rrit = 175GeV and mb = 4.7GeV. We have used 



the two- loop running MS quark masses |T2| and strong coupling constant ||13[, the CTEQ5M 
PDF and we choose the factorization and renormalization scale as m/^± + mw, if not 
stated otherwise. The expressions given in the previous section are valid for a general Two- 
Higgs-doublet model; they cover the MSSM case for a specific choice of the Higgs-boson 
masses and mixing angles. Here we focus on the MSSM scenario, taking into account the 
one-loop relations [|l^ between the Higgs-boson masses MhH,A,H^ and the angles a and 
/3, with 171^+ and (3 chosen as the two independent Higgs-input parameters, together with 
Ms = 1 TeV as a genuine SUSY mass scale. As a remark, the SUSY-QCD corrections 
arising from virtual gluino and squarks may also become important for specific parameters, 
which will be studied separately. 



In Fig. 2 we present the total cross sections, for the 66-annihilation subprocess, in the MS 
and OS renormalization schemes, and the relative deviation A = {cx^'^ — a'^^^) / {a^'^ + cr^^*^) 
in Fig. 3, both as functions of mjj±. The results differ remarkably at tree-level (LO) and 
are significantly closer at the level including the QCD corrections (NLO). Whereas in the 
OS scheme the LO cross section is sizeably reduced at NLO, the cross section increases in 



the MS scheme at NLO. 

For tan/5 = 50, A is still relatively large, around 8%, which can, however, be easily 
understood as a higher-order effect. For large tan/3 we can write the cross section in the 
two schemes as 

a0^ocm2p,,,(l-A + i?), 

a^ocm,(/i)2(l + 5), (19) 

with 



^ = 2^(log4 + t)' (20) 
71 mi 'i 
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which is the quantity entering also the relation between the MS quark mass and the corre- 
sponding pole mass at one-loop order, 

^(/^f =^?o.e[l-2-(log4 + |)]- (21) 

Hence, one finds 

A~^S/2. (22) 

For /X = 200GeV", we have A ~ 0.6, B ~ 0.3 and A ~ 9%. In the following we will give 
the results only in the MS scheme, which is usually more stable with respect to missing 
higher-order terms. 

Fig. 4 and Fig. 5 show the LO tree-level cross section and the NLO cross section 
(jNLO ^ including the hg contribution, as well as the relative correction 8 = [a^^^ — (T°)/(T° 
as a function of mH+ and tan /3. As pointed out before, due to the competition between 
the top-Higgs and bottom-Higgs Yukawa couplings, the cross sections are relatively small 
for intermediate values of tan/3, around tan/3 ~ 6 (see Fig. 6). From the figures, we can 
see that the QCD corrections are negative and the magnitude is greater than 14% for all 
charged Higgs mass and tan /I For tan/j = 2, the QCD corrections can decrease the cross 
section more than 30%. Table 1 contains numerical results for the relative correction 5 for 
a low and a high value of tan (3, tan (3 — 2 and 50. 

Fig. 6 and 7 display the LO and NLO hadronic cross section and the relative correction 
5 as a function of tan /3, for mH± — 200, 500 and 1000 GeV. From the figures, one can see 
that the QCD corrections are almost independent of tan/j when tan/3 > 15. For tan/3 < 15, 
the magnitude of the QCD corrections decreases with the increment of tan j3. 

In Fig. 8 we present the uncertainty of varying the factorization scale, which is defined 

as 

a(/x/ = 2/xo) - = |//o) 

with /Iq — rriHi + inw, as a function of m//±. We can see that the uncertainty is less than 
10% for all the charged Higgs boson mass and tan /3. 
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IV. CONCLUSION 



In conclusion, we have calculated the 0{as) corrections to the cross sections for bb 
W^H^ associated production at the LHC in both MS and OS renormahzation schemes, 
with an explicit discussion for the charged Higgs boson of the MSSM. The results in these two 
schemes are in good agreement. Our numerical NLO results include also the contributios 
from gb — > bH'^W^. In the MS scheme, the QCD corrections are negative and within 
-14% ~ -20% for charged Higgs mass up to 1 TeV and tan/3 > 15. For tan/3 = 2, the 
magnitude of the QCD corrections can be greater than 30%. 

The analytical results given in this paper are also vahd in a general 2-Higgs-doublet 
model, where the constrains among Higgs masses and angles a and f3 are released. As a 
final remark, the higher-order contributions arising from the gg — > H^W^bb, may also be 
important, especially for lower values of the charged Higgs-boson mass. They are presently 
under investigation. 
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VI. APPENDIX 



In this Appendix, we will present the non-vanishing form factors in the decomposition (|^ 
of the various s- and t-channel contributions to the one-loop matrix element. The form 
factors contributing to 6M^^^^ are 

4:^/2^T{mf — t)H 
+ (m^ - 6mli + P)Bo{i, 0, m^)], 

asghtUlt cos (3 2 , 2d n 2\ , / 2 , r\ D n 2m 

2v2vr(mt — t)^ 
/a = /iMfe, 

/g = (24) 
The form factors that contribute to 5M^^^^^ are (the argument the C functions is 



4v27r(m^ — t) 



ooj 



/2 = yJi'TTi ^-^ + 25o(t, A^ ml) + 4m^Co + 2(m,2 + ml - i)C, + 2m^C2 - 4Coo], 
4v27r(m^ — t) 

^3 = T^^T^^fl - 2i?o(t, A^ ml) - 2mlCo - 2{ml - i)C, - 2mlC, + 4Coo], 
A\/2Tx[mi — t) 

_ a sghtmtmb cos P 

2V2nim^-i) ^^'^^'^^ 

h = -^^iH^t-l + 25o(t, A^ m?) + 2mlCo - 2{ml + m? - m^)Ci + 2m^C2 

-4Coo-2m2Cn-2(m2-£)Ci2], 
_ asghtmtmbcosp ^^ 

V27r(m^ — t) 
_ a^^/ifembsin/? 
\'27r{mf — t) 
asghtmt cos pCi2 

/i2 = r- 2 — I^- (.25) 

V27r(m| - t) 

The form factors that contribute to 6M^°^ are (the argument of the C and D functions 

is [mlmls,ml,X^,ml] and [mj^±,ml,mlml,i, s,ml,ml, X'^,ml]) 

fi = ^""^^^ [htmbmt cos I3{Dq + D2) + hbsm.l3{CQ + m^Z^o + m'jj±Di 
2v 27r 
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+{mjj± +ml- u)D2 - {ml - u)Ds)], 
f2 = 2^J ^^^^ cos/3(-Co + irnl - ml - u)Dq + {ml - m^± - u)Di 
+ (m^ - m]j± - ml)D2 + {ml - m^ - u)D3 + 2Doo) 
+hbmbSmP{—Co — m^Do + {ml — u)Di — iD2 
-mlD^ + 2Doo)l 

h = ^^[/itmftmt cos/3(L'o + + ^2 + D^) + hsin(5{ml{D^ + D^) - {2ml - ^)^2)], 
/4 = -^^[htmt COS (3{Do + D, + D2 + D3) + hbmbsm(3{D, + D3)], 



/s = [htm^mt cos + /i^ sin (5{2Co + £)o + (Ci + C2) 



:7i 

-m^'Di + ^2 + mlD^ - 2L>oo + {ml - 2m^±)L>n 



h = -^^[/itmtCOS/3(L>o + 2L'i + £(2 + ^11 + ^12) + /i^mb sin /3(L>i + L>n + 1^12)], 
_ asghmbPs sin 

h = cos/5(Do + 2Di + ^2 + 1^3 + ^11 + ^12 + A3) 

V 27r 

+/ibm6sin/3(Di + + D12 + i^is)], 

/g = -^^[^tm;,mt cos /3(Dii + L>i2) + /i6 sin /3(2Co + m^(Do--Di) 
V 27r 

+C2 + (m^ - + - ^v^^^z + m^Dg - 2Doo 

+(m^ - 2mjj±)Dn + {2ml - 3^l± - ^1^12 

+{2ml -i- u)Di3 + {ml - mjj± - t)D22 

+ {2ml-i-u)D2^)l 

f^^ ^ _'^I^[htmtCOs(5{Do + 2^1 + 2^2 + Ai + 2r>i2 + ^^22) 
V 27r 



+ sm phbmb{Di + D2 + Dn + 2Di2 + D 



22) 



Ois9_ 

s/2'K 

hmb sin /3(L>2 + Dn + 2D^2 + i^ia + ^^22 + ^^23)]- (26) 



/12 = ^[^t^t COS/3(Di + ^2 + i^ll + 21^12 + 1^13 + i^22 + i^23 

V 27r 



The form factors contributing to 5M^'^'^^\A) are (the argument of the C functions is 
[ml,ml,s,ml,\^,ml]) 

h = h = ;r^^^^^^^[-l + ^B,{ml ml \') + 2mlC, + sC^ + 25C2], 
2v27r(s — m^o) 
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/6 = /lo = -/s. (27) 

The form factors which contribute to 6M^^'^'^\H) are (the argument of the C functions 
is [ml, ml, s, ml, X^, ml], the form factors to 5M^^^^\h) can be obtained by substituting 
cos a sm{a — P) — sin a cos(q; — /?)) 

75 = 76 = 79 = 710 = ^ /- ^- [1 - 2So(mfe, m^, A ) 



2^2' 



vr s — m 



-AmlCo - {2ml + - 2SC2]. (28) 

The form factors which contribute to (JM^^^*^ are (the argument of the C functions is 

[i, ml,ml±,mlX^,ml]) 

/i = 0(sgmb j /^^^^^^^og^(-(g^ _^ _^ /ifcSin/3(-l + 2Bo{ml,ml, A^) 
2v27r(t — ) 

+3mj2(7o + {-ml + mj^± + m^ + i)Ci + {-ml + 2m^± + m2)C2)], 
2v27r(t — mf) 

+ {ml - m^± - 2i)Ci + {ml - 2m\± -i)C2) + hmt sin /3{-mlCo 
-iC^ + {ml~i)C2)], 
/g = 2/3 = 2/i/mfe. (29) 



The scalar integrates B,C, D in above expressions are evaluated by FF package |16 . 
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FIG. 1. Feynmann diagrams for the subprocess bb ^ H at the NLO: (a)-(b) Born diagrams; 
(c)-(h) virtual correction diagrams; (i)-(m) gluon-radiation diagrams. 
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Charged HIggs Mass (GaV) 



FIG. 2. The total cross sections of process PP — > 66 — *■ H^W^ versus mij± at the LHC with 
^/s = 14 TeV in MS and OS schemes. The hnes correspond to the cross sections at the one-loop 
level of MS (solid), one- loop level of OS (dashed), tree level of MS (dotted) and tree level of OS 
(dash-dotted). The three groups correspond to tan/3 = 50, 2, 6 from above to the bottom. 
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FIG. 3. The A (defined in text) versus mfj±. The thick hnes stand for the results of one-loop 
and thin lines stand for the results of tree level, tan/? = 2 (solid), 6 (dashed), 50 (dotted). 



17 




Charged HIggs Mass (GaV) 



FIG. 4. The total cross sections for the subprocesses bb — *■ H^W^ and gb(b) H^W^b{b) 
versus m^± at the LHC with ^/s = 14 TeV. The thick hnes correspond to the one-loop cross 
sections with tan/3 = 2 (solid), 6 (dashed) and 50 (dotted). The tree level cross sections are also 
shown by the corresponding thin lines. 
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FIG. 5. The QCD relative corrections to the cross sections for the subprocesses bb — *■ H^W^ 
and gb(b) ^ H'^W^bib) versus mH±- The Hnes correspond to tan/? = 2(sohd), 6 (dashed) and 
50(dotted). 
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FIG. 8. The uncertainty [defined in Eq. (|23|) ] as a function of mjj±, where the renormahzation 
scale is fixed as mjj± + rriY/. The hnes correspond to tan/3 = 2(solid), 6(dashed) and 50(dotted). 
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mH± (GeV) 


tan 13 = 50 


tan ^ = 2 


200 


~ -16% 


- -32 % 


500 


~ -15% 


~ -34 % 


1000 


~ -14% 


~ -30 % 



TABLE 1. QCD corrections to the process bb H^W^ in MSSM and 2HDM in MS scheme. 



The renormahzation and factorization scales are taken as mfj± + mw- 
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